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Abstract. Microquasars are binary systems that harbor a normal star 
I and a compact object (black-hole or neutron star), and show relativistic 

outflows (or jets). The matter that forms these jets is of likely stellar 
origin, previously expelled from the star and trapped in the potential well 
of the compact object. This matter is accreted by the compact object, 
forming a disk due to its angular momentum, and is eventually ejected in 
' the form of a bipolar outflow (the jets), which generates radio emission 

and could also be a very high-energy emitter. To study and understand 
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O ' the radiation from microquasars, there is a set of elements that can play a 

+l3 [ major role and are to be taken into account: the photons and the expelled 

^ ' matter from the star in the case of high-mass systems; the accreted matter 

radiation; the jet; the magnetic field carried by the jet or filling the bi- 
nary system; and the medium surrounding the microquasar at large scales 
(~ pc). In this lecture, we consider these elements of the microquasar sce- 
1^ ' nario and briefly describe the physical conditions and processes involved in 

. the production of non-thermal radiation from radio to gamma-rays. The 

' required energetics, particle acceleration and transport, several radiative 

mechanisms, and the impact of different photon absorption processes, are 
' discussed. 

o 
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Q , 1. Introduction 

Among the different classes of astrophysical sources, microquasars (see Mirabel 
^X^ ■ k Rodri guez 1999 for a review) are specially interesting for the study of different 

topics of high energy astrophysics. The fact that these objects present non- 
thermal emission in different wavelengths, from radio (e.g. Ribo 2005) to infrared 
wavelengths (e.g. Mirabel et al. 1998), and from X-rays (e.g. Corbel et al. 
2002) to very high-energy (VHE) gamma-rays (Aharonian et al. 2005, 2006a; 
Albert et al. 2006; Albert et al. 2007), has several implications, being perhaps 
the most important the following ones: particle acceleration takes place, which 
means that there are large amounts of energy contained in moving matter and 
magnetic fields in a low entropic state, ready to be released in the form of heat and 
radiation; the acceleration efficiencies are high, pointing to quite extreme plasma 
conditions from the point of view of particle confinement, flow velocities, and 
specific characteristics of radiation, matter, and magnetic fields in the emitting 
region; the non-thermal energy outcome can be very large, which implies that 
the main radiative mechanism is efficient (e.g. close to the saturation regime, 
when practically all the particle energy is radiated away). In addition, extended 
emission in the radio (e.g. Mirabel et al. 1992) and X-ray bands (e.g. Corbel 
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et al. 2002) has been detected, in some cases up to very large scales, implying 
that there is still efficient non-thermal particle acceleration at the sites where 
this radiation is generated. 

Within the class of microquasars, the subclasses of low- and high-mass sys- 
tems are distinguished depending on the mass of the stellar companion, which 
can be smaller or bigger than the compact object mass. In general, high-mass 
microquasars have OB stars as companions, which produce strong stellar winds 
and are very bright in the optical/UV band. Therefore, unlike low-mass mi- 
croquasars, whose radiative processes would be determined basically by the ac- 
cretion/jet system only, high-mass microquasars should present a more complex 
phenomenology due to the presence of the wind and the radiation field from 
the star. In such an environment the non-thermal radiation would be affected 
in different ways, being perhaps the most relevant ones: dynamical interaction 
between the jet and the stellar wind (e.g. Perucho & Bosch- Ramon 2008); atten- 
uation of the radio emission due to free-free absorption in the wind (e.g. Szostek 
& Zdziarski 2007); inverse Compton (IC) scattering of relativistic electrons with 
stellar photons (e.g. Paredes et al. 2000; Kaufman Bernado et al. 2002); the 
stellar wind could provide with targets for proton-proton collisions (e.g. Romero 
et al. 2003; Aharonian et al. 2006b); the gamma-rays of energy above the pair 
creation threshold would be absorbed via photon-photon interactions with the 
radiation field of the star (Boettcher & Dermer 2005); and the created secondary 
pairs would radiate under the magnetic and radiation fields present in the star 
surroundings (Bosch- Ramon et al. 2008). 

In this lecture, we try to give a broad overview of the physical processes 
involved in the formation and characterization of the observed non-thermal emis- 
sion from microquasars. In Section [2?1 a qualitative and semi-quantitative de- 
scription of the physical scenario is given, and the high energy processes relevant 
in microquasars are introduced; in Sect. [3], the variability and spectral proper- 
ties of the emission are discussed; finally, in Sect. |4]1 we conclude summarizing 
plus some comments and remarks. 

2. High energy processes in microquasars 
2.1. The physical scenario 

In Fig. [H we present the picture of a microquasar including the elements of the 
scenario: the star, the accretion disli^, the compact object, the jet, the magnetic 
field, a relativistic population of electrons and protons in the jet, the stellar 
wind, the main radiative processes, and pair creation. Rather than going to 
details about the processes of generation of jets, magnetic fields, acceleration 
of particles, and the like, we will give for granted what we already know from 
observations, i.e. that microquasars present all these ingredients. Instead of 
that, what can be done is to constrain the models that describe these physical 



^The accreted matter has angular momentum, which leads to the formation of a disk. Some 
mechanism (e.g. turbulent viscosity -Shakura & Sunyaev 1973- or an outflow -Bogovalov & 
Kelner 2005-) can remove energy and angular momentum, allowing the material to drift towards 
regions of the accretion disk closer to the compact object. 
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entities or processes using observations and basic tools from elementary physics. 
These tools are presented in the next sections. 

2.2. Energetics 

Non-thermal radio emission is produced in the jets of microquasars, but the 
energies of the photons produced in these jets likely go up to gamma-rays (e.g. 
Paredes et al. 2000). Therefore, the first thing to do, in order to understand the 
microquasar non-thermal emission, is to look for the origin of the energy powering 
these jets. Natural energy sources are the potential energy of the accreted matter 
(see sect. 3.1.2 in Bosch-Ramon et al. 2006 for a semi-quantitative discussion), 
or the rotational energy of the compact object (e.g. Semenov et al. 2004; see 
also Li et al. 2008). Since jet activity seems to be well correlated with accretion 
in microquasars (Fender et al. 2004), we will focus here on the accretion energy 
requirements. The observed non-thermal luminosities can reach large values, 
~ 10^^—10^^ erg s~^, and the jet kinetic luminosity is, order of magnitude, similar 
or above these values, i.e. ^ 10^^ erg s~^. For an efficiency of the accretion- 
ejection energy transfer of a 10%, the required accretion energy budget will be 
^acc ~ 10'^'' erg s^^. This value is quite modest, being a few % of the Eddington 
luminosity for a few M© compact object. In the phenomenological classification 
of Fender et al. (2004), this Lacc-value would correspond to a microquasar in the 
low-hard state, when a persistent jet is expected to be present, i.e. microquasars 
can persistently power gamma-ray emission. 

2.3. Particle acceleration 

Presently, it is very difficult to distinguish between several mechanisms of par- 
ticle acceleration in microquasar jets (for a discussion of this, see Bosch-Ramon 
2007 and references therein). Due to this fact, for the purpose of this lecture, 
it is more convenient to investigate the acceleration of particles using a simple 
approach. The first step is to see whether a necessary condition for the accel- 
eration of a charged particle with certain energy is fulfilled, which is given by 
the Hillas criterium (Hillas 1984). This consists on the fact that particles can 
only be accelerated if their Larmor radius (rL = E/qB^; where is the accel- 
erator magnetic field, and q and E are the charge and energy of the particle, 
respectively) is smaller than the accelerator size (la); otherwise particles escape 
the accelerator. This limits the highest achievable energy to (if not specifically 
stated, the units in this work are cgs): 

E < qBJ., ■ (1) 

Making a step further, to determine whether particles can be accelerated up to a 
certain energy, the specific acceleration and energy loss (or particle escape) mech- 
anisms are to be known: t^cc = ^cool/esc • general, the acceleration timescale 
can be expressed as: 

iacc = ri— = r]—— , (2) 
c qBc 

where rj is a dimensionless phenomenological parameter (or function) represent- 
ing the acceleration efficiency, different in each acceleration scenario. The par- 
ticular case of 77 = 1 corresponds to the shortest possible acceleration time in- 
dependently of the acceleration mechanism. An instance for 77 can be given for 
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Figure 1. Illustrative picture of the microquasar scenario, in which 
the main elements and processes considered in this work are shown 
(background image adapted from ESA, NASA, and Felix Mirabel - 
CEA, IAFE/CONICET-). 
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the case of non-relativistic diffusive shock acceleration (plane shock with weak 
magnetic field, in the test particle approximation -Drury 1983-): 

where V^h is the shock velocity, and D is the diffusion coefficient (-DBohm in the 
Bohm limit). For = 3 X 10^ cm s"^ and D = L>Bohm, r/ ~ 10^. 

Without focusing on any particular acceleration mechanism, it is worthy 
noting that the final product of particle acceleration of a diffusive or stodiastic 
mechanism will be a power-law {Q{E) oc E~^), and /i ~ 2 is commonly adopted. 
Diffusive non-relativistic strong shock acceleration results in /i = 2. Another 
point to notice is that, usually, the particles in the emitter, e.g. the jet, are 
considered to have an isotropic velocity distribution in the jet reference frame 
(RF), due to deflection in tlie randomic component of the magnetic field, which 
moves solidary with the jet matter. 

2.4. Particle propagation 

In a medium of diffusion coefficient D, particles propagate in one particular 
direction witli the time dependence Idm ^ ^/2Dt. Under tlie impact of cooling, 
the typical distance that particles of energy -ExeV = E/1 TeV, and under a 
magnetic field of Bq = B/1 G, can reach is: 

W-10-<^i?-/^*;/2l(7y^)'''cm, (4) 

V-t^Bohm/ 

where the Bohm rate is the slowest possible diffusion rate, i^Bohm = '"lc/S, and 
^cool is the cooling timescale in seconds of the dominant loss mechanism. If the 
medium in which particles are embedded is also moving (e.g. a jet), there is in 
addition advective transport: 

These formulae allow us to estimate, under different cooling processes, whether 
particles can propagate significantly wlien comparing these distances with the 
accelerator size. For instance, jet advection may transport particles up to large 
distances. If the radiative processes were still efficient, such a situation would 
imply that the emitter size Z 3> Za (for a thorough discussion, see Khangulyan et 
al. 2008). 

2.5. Cooling processes 

We consider here the following non-thermal processes: relativistic Bremsstrahlung, 
synchrotron and IC (for an exhaustive review, see Blumenthal &; Gould 1970), 
when leptons are the main emitters; and proton-proton interactions (see Kelner 
et al. 2006), wiien protons are the particles that produce the emission to study. 
Here we briefly list the relevant cooling timescales of these processes. The adia- 
batic loss time is also given, since relativistic particles could lose energy exerting 
work on the medium. 
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Relativistic electrons radiate Bremsstrahlung under the effect of the electric 
field in the surroundings of an atom. The characteristic cooling timescale is (e.g. 
Cheng & Romero 2004): 

^Br ' S , (6) 

n 

where n is the density of the medium. Synchrotron radiation originates when 
an electron moves in an irregular magnetic field, spiraling around the chaotically 
oriented magnetic lines. The electron suffers Lorentz forces and radiate in the 
direction of motion the energy corresponding to the momentum perpendicular 
to the magnetic field. The typical frequency of the outgoing synchrotron photon 
is V = 6.27 X lO^^BE^. It is thought that a significant part of the magnetic 
energy density in jets is in the form of a chaotic magnetic field (e.g. Biermann 
& Strittmater 1987) attached to the plasma. This random magnetic field also 
confines the relativistic particles, which diffuse inside the jet instead of escaping 
from it almost at the speed of light. The characteristic energy loss timescale for 
synchrotron emission is: 

tsy ~ 4 X W^B^^E^^y s . (7) 

Under the presence of radiation fields, relativistic electrons also interact with the 
ambient photons via IC scattering, in which the scattered photon has increased 
its energy by 4(£'/meC^)^ and takes the direction of the scattering electron. If the 
energy of the target photon (eo), in the electron RF, is < mgC^, the interaction 
takes place in the classical Thomson regime; in the emitter RF, the electron 
energy must be < mgC^/eo. The associated characteristic timescale is: 

16 

^IC T 



Wrad-^TeV 



where Urad is the target photon field energy density (cgs). In case the target 
photon, in the electron RF, has an energy > rrieC^, the interaction enters in the 
Klein-Nishina regime, with a characteristic timescale (Khangulyan et al. 2008): 

tic KN ~ lO^nrad)-' (j^) ^ ^ ^TeV S ■ (9) 

It is worthy noting here that IC scattering has a strong dependence on the angle 
of interaction between the incoming electron which becomes less important when 
entering in the Klein-Nishina regime. This has to be accounted when computing 
the radiation for a specific system and orbital phase, since the geometry of the 
interaction changes along the orbit (see Khangulyan et al. 2008 for a deeper 
discussion on this; see Blumenthal & Gould 1970 for the angle-averaged cross 
section, and Bogovalov & Aharonian 200CE] for the angular dependent one). We 
note however that the particle distribution is usually assumed isotropic in the 
emitter RF and tic will not depend on the angle. This emitter may be moving 



^Eq. (19) in this work lacks a (1 — cos(6')) inside the integral {6 is the angle between the electron 
and photon directions of motion in the laboratory RF). 
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in a certain direction with respect to the observer, e.g. in the jet, in which case 
Doppler boosting has to be taken into account. 

If relativistic protons are present, they will lose energy interacting with the 
ambient atoms (e.g. Romero et al. 2003; Aharonian et al. 2006b). At the 
relevant energies, i.e. > GeV, ionization losses become negligible, and proton- 
proton collisions are the dominant cooling channel with timescales: 

10^5 cm-=^ 

*pp~ (10) 

nt 

where rit is the density of targets (cgs). After colliding, among other possible 
interaction products, about half of the proton energy goes to 7r°/^, which decay 
to gamma-rays///^ /f^, and /U^ decay to (see Orellana et al. 2007 for their ra- 
diative relevance in microquasars), Uc and f^. Roughly, 1/6 of the proton energy 
goes to gamma-rays, a similar amount to u, and 1/12 to e^. If dense and very 
hot photon fields and ultrarelativistic nuclei were present, photo-meson produc- 
tion and photo-disintegration could play some role, and even synchrotron proton 
radiation may be efficient under certain conditions (Vila & Romero 2008). Here 
we will not consider them for simplicity, but also because they are in general little 
efficient or start to work at energies that may not be reached due to acceleration 
constraints (Bosch-Ramon & Khangulyan 2008, in preparation). 

The cooling timescale for particles exerting work on the surrounding medium 

is: 

21 , , 

iad~:TT^s. (11) 

"J ''cxp 

This can be considered to happen in case the relativistic particles are confined 
within a medium that suffers expansion, like a jet in overpressure with its envi- 
ronment. 

2.6. Photon attenuation processes 

Absorption of gamma-rays with energies above the pair creation threshold (eth = 
2mgC^/eo(l — cos(^?)); where is the target photon energy and Q is the interaction 
angle between the directions of the VHE and the stellar photons) will occur if 
there is a photon field surrounding the VHE emitter. In case of microquasars 
harboring massive stars, the opacity for one photon to be absorbed, related to 
the interaction probability, can be close to or much larger than 1, i.e. photon- 
photon absorption is either optically thin or optically thick. To estimate the 
significance of this process, we can use the next simple expression: 

^^^^ 10-25 (12) 

where typically the system size i?oib ~ 10^^ — 10^^ cm, ~ 1.6 x 10"^^ erg (T* k, 
4 X 10^ K) is the typical energy of the stellar photons, and tirad ~ 100 erg cm"^. 

is the opacity for gamma-rays with energies ~ 3.7 x cth in an isotropic target 
photon field (when the cross section is the largest; see Coppi & Blandford 1990, 
eq. 4.7). 

Actually, in most high-mass microquasars (LS 5039, Khangulyan et al. 2008; 
Cygnus X-1, Bednarek & Giovanelli 2007; Cygnus X-3, Protheroe & Stanev 1987; 
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SS 433, Reynoso et al. 2008; and LS I +61 303, Romero et al. 2007) opacities 
significantly above 1 will occur. The occurrence of absorption is important not 
only because it reduces the amount of VHE photons that escape the system, but 
also because it depends on 6. This dependence is present through the interaction 
probability and eth- In microquasars, since the interaction angle changes along 
the orbit, as well as the surface density of target photons seen by the observer, 
opacity changes along the orbit (a discussion of this can be found in Khangulyan 
et al. 2008; the angle averaged pair creation rate can be found in Coppi & 
Blandford 1990, and the angle dependent cross section in Gould & Schreder 
1967)0. 

Electromagnetic cascade can take place efficiently when KN IC losses are 
dominant and r^^ ^> 1. The former implies that the ambient magnetic field 
outside the emitter must be low enough, i.e.: 

^erg s" 



< 10 r — G. (13) 

\ 1039 erg s-i; \rJ ^ ^ 



At i?c, KN IC losses are equal to synchrotron losses for 1 TeV photons. Electro- 
magnetic cascading consist on a VHE electron that creates a VHE photon via 
IC, which is absorbed in the stellar field creating a VHE electron/positron paii0, 
which to its turn creates another VHE photon via IC, etcetera (this is quite 
constraining, actually, since could be significantly larger than Be, as noted in 
Sect. 12.71) . Electromagnetic cascading occurring deep in the KN regime reduces 
effectively the opacity of the system. If the magnetic field were similar to or 
larger than B^ the absorbed energy would be mainly released via synchrotron 
radiation (see Bosch-Ramon et al. 2008). 

Beside gamma-ray absorption due to pair creation, it is worthy also to note 
that radio emission can be attenuated due to free-free absorption in the stellar 
wind. An estimate for the free-free opacity (Rybicki & Lightman 1979) can be 
obtained from the next formula: 

ff ^ 40- , (14) 



7y2 P3 T/2 rp. 

'^5 GHz-f^O.3 AU'^w S.S-'^ 



3/2 
4 



where M7 = M/IO"'^ Mq yr ^ is the stellar mass loss rate, 1^5 gHz = v/^ GHz the 
frequency, i?o.3 = -R/0.3 AU the distance to the star, 8.3 = Vw/2 x 10® cm s~^ 
the stellar wind velocity, and 4 = T^/IO^ K the wind temperature. This 
estimate could be affected by the real structure of the stellar wind, which may 
be clumpy (e.g. Owocki & Cohen 2006), reduced or increased depending on 
the clump mass density and number. This estimate suggests that the observed 
radio emission is produced outside the binary system in high-mass microquasars 
(^ AU), at spatial scales that may be only marginally resolvable by the present 
radio interferometer instruments, but still far from the jet ejection region. 



■^When computing the opacity, the factor (1— cos(0)) should be also considered in the calculation, 
as in the case of angle-dependent IC scattering. 

*Well above the threshold, when cascading is efficient -deep KN IC-, one member of the pair 
will take most of the initial photon energy (e.g. Boettcher & Schlickeiser 1997). 
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Another source of radiation attenuation is synchrotron self-absorption. Elec- 
trons of certain energies, within a source of certain size and magnetic field, can 
efficiently absorb the synchrotron emission produced by them. For an homoge- 
neous emitter with a population of relativistic electrons producing synchrotron 
emission, the synchrotron self-absorption frequency is (see chapter 3 in Pachol- 
czyk 1970): 

where K = N(E)/E~p [p here is the one at the relevant electron energies) is 
the electron energy distribution normalizatioE0, ci = 6.3 x 10^^, C2 ~ 10"^'' (c2 
actually depends on p, but the given value works for p f« 2 — 2.5), and I is the 
size of the emitter. We notice that, if the emitter has a complex structure or is 
strongly inhomogeneous, fgsa may be significantly different. 

Finally, ionization of the stellar wind atoms can absorb the lower energy 
part of the X-ray spectrum, and also induce relativistic particles cooling. We do 
not discuss these processes here though note that they may be relevant in some 
specific situations (e.g. Bosch-Ramon et al. 2007; 2008). 

2.7. Targets for radiation and energy losses 

Relativistic particles interact with different targets: magnetic, photon and mat- 
ter fields, producing emission and subsequently losing energy. In microquasars, 
there are different sources for the mentioned targets. Here, we will focus on 
those related to the star (for the case of high-mass systems), the jet, and the 
environment at large scales. In the regions close to the compact object, where 
the jet is thought to be produced, dense target fields from the accretion disk 
could be present as well, but given our lack of knowledge concerning the real jet 
and environment properties in those regions, we do not consider them here (for 
a brief discussion, see Bosch-Ramon 2007). 

The star is a very important ingredient in high-mass microquasars. On one 
hand, the stellar photon field, with typical luminosities ~ 10^^ — 10^^ erg s~^ 
and temperatures ~ 4 x 10^ K, renders stellar IC a very efficient process, as 
well as photon-photon absorption. On the other hand, the stellar wind provides 
targets for proton-proton collisions, could be a source of radio and soft X-ray 
radiation attenuation, and may play a dominant role on the jet dynamics at 
spatial scales ~ Rorh (not treated here; see Perucho & Bosch-Ramon 2008 for a 
thorough discussion). The wind can have a very complex structure as noted in 
Sect. 12.6.1 although we estimate its density assuming that this is an homogeneous 
supersonic radial flow of velocity ~ 2 x 10® cm s~^: 

o / \ /3 X 10^2 \ 2 



^It is worthy noting that the injection electron spectrum Q{E) and the particle energy distri- 
bution N{E) are different functions. The first one gives the energy dependence of the injected 
particles, and the second one gives the same but for the final particles that suffered the different 
cooling and escape processes that take place in the emitter. The units may not coincide. 
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Concerning the stellar photon field, the photon number density (valid for a black 
body) can be roughly approximated by: 



For the magnetic field in the surroundings of the star, we can adopt a R- 
dependence of -Ba in the range r = 1 — 3 (Usov et al. 1992), which would cor- 
respond to a toroidal-radial-dipolar dependence normalized to -6=,, at the stellar 
surface (-R*): 

It is worthy noting that is typically 100-1000 G in OB stars, found for instance 
using the Zeeman effect (Donati et al. 2002). Further evidences of the presence 
of the magnetic field come from the detection of non-thermal radio emission 
from isolated massive stars (e.g. Benaglia 2005 and references there in), and X- 
ray observations (e.g. Stelzer et al. 2005). Recalling what was said in Sect. 12.61 
regarding electromagnetic cascades, the typical -B,,-values would seem to preclude 
the occurrence of efficient cascading in high-mass microquasars. 

The magnetic field in the jet could be assumed to go like 1/z (e.g. for a 
toroidal magnetic field), where z is the jet height, and normalized to the magnetic 
field value at some height (i.e. Bq at zq): 



f Bo \ f 10^ cm\ / 10 



12 



cm 



The presence of magnetic field and relativistic electrons in the jet leads to syn- 
chrotron emission. This process generates a target photon field inside the jet 
that could be suitable, under certain conditions of jet radius (-Rjet), Rjet, and 
injected particle luminosity (i^rei), for efficient synchrotron self-Compton (SSC), 
i.e. scattering of the synchrotron photons by the relativistic electrons that pro- 
duced them. A relationship that tells when synchrotron self-Compton losses 
start to dominate over synchrotron ones is the following: 



^ict ^sync , ^rel 
""B = ^ —52— ' 20 



where the main assumption is that the synchrotron emission (of luminosity Lgync) 
comes mostly from one compact region of magnetic field -Bjet- The treatment of 
SSC losses, when dominant, is very complicated since it is non-linear, i.e. the 
radiation and the particle distribution are strongly coupled, with feedback effects. 
Therefore, it is important to know when SSC losses are to be accounted for (e.g. 
Bosch-Ramon & Paredes 2004). Actually, in some cases the SSC radiation can 
be computed without the need of taking into account its losses since other energy 
losses, like synchrotron itself, or external IC, are dominant. Finally, the matter 
density of a conical supersonic jet can be expressed in the next form: 

' Uet \ ( 10^° cm s'i\ / lO'^ cm^ ^ 



15 I jet \ / -LU ciii s \ / iu ciii \ _3 



njet = 2 X 10^^ I ^ :r^^^ r 1 I I I I Cm 



(Fjet- 1)1037 ergs-1 / \ V;,t \ R 



jet 



(21) 
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where Vjet is the jet speed and Fjet the jet Lorentz factor. In general, densities 
are expected to be low, probably not enough to produce gamma-ray emission via 
relativistic Bremsstrahlung, or proton-proton interactions (regarding the latter, 
see otherwise the case of SS 433, which has very heavy jets -Reynoso et al. 2008-). 

At the largest scales, where the microquasar jet terminates, the treatment of 
the ambient radiation, matter and magnetic fields becomes more complex. For in- 
stance, in massive systems, a powerful stellar wind would reduce the density and 
increase the pressure of the jet environment via shocking the ISM, modifying the 
whole dynamics of the jet/medium interaction and introducing inhomogeneities 
in the surroundings. Also, high-mass microquasars, being young objects, tend to 
be embedded in a dense medium, unlike low-mass microquasars, which can be 
quite old and present even in the halo. In the case of low-mass systems, the envi- 
ronment of the microquasar may have very different densities, 10~^-> 1 cm~^, 
depending on the location within the Galaxy, and the ISM, which could present 
by itself strong inhomogeneities at pc scales. In addition, a massive and hot star 
could generate a radiation field above the local one (i.e. the galactic background). 
Otherwise, the latter would be the only possible target for IC interactions, which 
would typically yield quite low efficiencies. Like the density, the magnetic field in 
the jet termination regions is an uncertain parameter, but some equipartition ar- 
guments and simple jet/medium interaction models give values ~ 10-5-10-3 G, 
as found e.g. by Bordas et al. 2008. In that work, the authors show that the 
non-thermal emission from the jet termination regions could be detectabl^. 

3. On the radiative output 

In this section, we present some examples of the spectra and the lightcurves that 
the emission from microquasars could show. In order to shape the spectrum 
and the lightcurve, several possible processes to explain observations at different 
wavelengths are considered. First, a semi-qualitative/quantitative insight into 
the formation of the particle energy distribution is given, as well as a brief de- 
scription of the shape of the spectrum of the radiation produced from different 
particle energy distributions. 

3.1. The energy distribution of particles and the radiation spectrum 

In Section 12.3.1 we have already mentioned that the injected particle spectrum 
is usually adopted with index h = 2. Although the complexity of the emitter 
(radiative losses, non radiative losses, escape losses, etc.) may require the use 
of a very complex model setup, the basic features of the evolved particle energy 
distribution can be extracted from a very simple case: a (one-zone) homogeneous 
emitter with particle escape in the steady regime (see Eq. 3-7 in Khangulyan et 
al. 2007). The main idea is that the final particle energy distribution, iV(£'jl|, 



Interestingly, microquasars are not the only galactic sources that produce non-thermal emission 
via jet-environment interactions (see Araudo et al. 2007 and references therein). 

^The emitter is considered steady and homogeneous, without spatial nor time dependences in 
N{E). 
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is: 

(X = E~^^^ / E > minimum injection energy , oc otherwise. (22) 

From E = — -E/tcooi E^ , we can roughly give the N{E) shape for different pro- 
cesses, where g = 2 in the case of synchrotron and Thomson IC, g' ~ in the case 
of KN IC and ionization losses, g ~ 1 in the case of relativistic Bremsstrahlung, 
g ~ 1 in the case of proton-proton interactions, and = 1 in the case of adi- 
abatic losses. If an escape time is introduced (e.g. due to convection by the 
outflow away from the shock that accelerates particles), no particles older than 
this escape time will be present, and they will not contribute to the particle 
energy distribution below the energy at which they escape. The particles that 
do not have time to cool will escape the source at the same rate whatever their 
energy. It implies that, in general, below the escape energy, N{E) oc E~^ (above 
the minimum injection energy) and = (below the minimum injection energy). 
Analogously, in the time dependent case, i.e the age of the accelerator < the 
cooling timescale at a certain energy, the particle energy distribution has p = h 
below this energy (so-called break energy). All this is valid for protons and for 
electrons. In Fig. [2l the energy distribution of electrons evolving under syn- 
chrotron cooling is shown to illustrate the effect of energy losses in the injected 
particle spectrum. 

For the shape of the spectral energy distribution (z^F(z^), where F(i') is the 
specific flux), it can be written: 

vF{v) = n{v) = v N{E)P{E, v)dE « vE^^{E{v))N{E{v))dE{v) jdv , 

(23) 

in the delta function approximation for the specific power -P(-E, v) of one par- 
ticle of energy -B, approximately valid for the radiative processes presented in 
Sect. 12.5.1 n{y) is the specific photon number. The cooling and the radiation 
rates, -Ecooi and -Eradj refer to the radiation we calculate and the dominant cool- 
ing mechanism affecting particles, respectively. E^^A and -Ecooi do not necessarily 
have to be the same. 

To go further, E(y) is required, and a dependence of the kind E(y) oc will 
be adopted, which is roughly correct for the mechanisms considered here. For 
synchrotron and Thomson IC, / = 1/2; for KN IC, relativistic Bremsstrahlung 
and proton proton interactions, / ~ 1. Therefore, for a cooled particle population, 

vF{v) oc vE,^{E{v))E{vf-''v'-^lE,,,,{E{v)) oc z.'(^^^ad-9cooi+2-h) ^ (24) 

which implies, for instance, that vF(y^ oc constant for synchrotron, Thomson 
IC and proton proton radiation under E'rad = -Ecool, or v~'^ for KN IC under 
synchrotron cooling. This second case is presented in Fig. [3l For an uncooled 
particle population, 

vF(y) oc vE,^{E(y))E(y)-^v^-^ oc i/Hff.ad+i-/*) ^ (25) 

and below the minimum injection energy (if particles had time to cool before 
escaping or at the age of the source) : 



vF(v) a vE,^{E{v)y-^lE,,,x oc z.'(9-d-9cooi+i) 



(26) 
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Figure 2. The sketch of a particle energy distribution times E'^ is 
shown. The electron population suffers synchrotron losses and the 
cooled particle distribution is oc X E-^ for h=2. In the case shown 
here, below the break energy the age of the source equals or is smaller 
than the synchrotron energy loss timescale, and the spectrum is still 
uncooled. 
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h=2; Thomson/KN IC SED under synchrotron cooling 
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Figure 3. The sketch of a spectral energy distribution of IC emission 
produced by electrons cooled by synchrotron emission. The synchrotron 
component is not shown. 
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Below the escape energy and the minimum injection energy, the delta func- 
tion approximation does not work, and the proper emissivity function is to be 
accounted for (e.g. for synchrotron /Thomson IC emission: vF{v) oc z/^/"^; for 
relativistic Bremsstrahlung: vF{v) cc v^). Complex shapes for the particle spec- 
trum, very different from a power-law, and the presence of low energy and high 
energy cut-offs or extreme slopes would modify the simple schema presented here. 



3.2. Variability and dominant emission mechanisms at difTerent wave- 
lengths 

In this section, we describe briefly the variabihty of the source depending on the 
energy range. For simplicity, we keep our discussion in the context of one-zone 
or homogeneous models. We note that propagation effects, which can lead to 
emitter inhomogeneity, can change the variability pattern of the radiation, e.g. 
smoothing it (see Fig. 22 in Khangulyan et al. 2008 for an example of the impact 
of particle advection in the lightcurve and spectrum), or increasing it via particle 
escape (e.g. Khangulyan et al. 2007). 

In the simple context adopted here, e.g. one-zone model, different factors 
can influence on the emission behavior, like the particle injection rate, the mag- 
netic field, the target photon density, the non-radiative cooling rate, the escape 
timescale, and the geometry changes of the binary /observer system along the 
orbit. It is worthy noting that in the saturation regime, the change of target 
densities does not necessarily render variabihty as long as the injected lumi- 
nosity is constant (e.g. changing the stellar photon density along the orbit in an 
eccentric system does not imply modulation of the IC emission). 

Radio 

Radio emission, of synchrotron origin in microquasars, can be affected by 
all the mentioned factors that influence on the lightcurve but the system geom- 
etry. Nevertheless, it is worth to note that synchrotron self-absorption makes 
the emitter to become opaque to its own synchrotron radiation below a certain 
frequency. In such a regime of emission, variations in the observed spectral break 
frequency (fssa; the optically thick-optically thin transition) can occur depending 
from which direction the observer sees an (e.g. asimetric) emitter. We also note 
that radio photons are subject to free free absorption, very likely to occur in the 
stellar wind. This will also introduce orbital dependent variability. 

X-rays 

X-rays could be of synchrotron origin, produced by the high energy part 
of N{E), or IC origin, produced by the low energy part of N{E) (i.e. in the 
Thomson regime). In the former case, the radiation would be affected by all the 
factors mentioned above but changes of geometry, unless the electron distribution 
and/or the magnetic field configuration are themselves anisotropic. Synchrotron 
self-absorption is hardly relevant at these energies. IC can be affected, together 
with the other variability sources, by geometrical effects via the next dependence: 
vF{y) (X (1 — cos(^)), since {9) will change strongly along the orbit since the 
stellar photons come from the star direction, which is orbital phase dependent 
from the observer point of view. Relativistic Bremsstrahlung can be hardly 
efhcient in this energy range, since it is very hard below the energies around the 
electron rest mass one (~ MeV) and will fall probably below the synchrotron 
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and/or IC fluxes. X-ray attenuation due to ionization of wind atoms, with a 
changing hydrogen column density along the orbit, may be a source of variability. 
Gamma-rays (GeV) 

The emission in this energy range could be due to relativistic Bremsstrahlung, 
proton-proton interactions, or IC scattering. If hadronic and the leptonic accel- 
eration were equally efficient, and IC scattering negligible (e.g. under a low-mass 
companion photon field) , relativistic Bremsstrahlung would dominate due to the 
smaller energy transferring in the case of (proton-proton) proton to gamma-rays 
compared with (relativistic Bremsstrahlung) electron to gamma-rays. Neverthe- 
less, very high densities would be required (see Sect. 12.51 and 12. 7T]) . If a massive 
and hot primary is present, stellar (Thomson) IC will be a much more efficient 
process. SSC radiation could be still relevant for the case of a high-mass micro- 
quasar, and certainly is a good candidate for low-mass microquasar^. As noted 
above, powerful jets containing relativistic protons in low-mass microquasars 
may also present significant hadronic radiation from other channels. 

The variability in this energy range could be due to all the effects mentioned 
above. In particular, the angular dependence of (Thomson) IC (~ (1 — cos(^))), 
implies that less IC emission goes to the observer when the emitter is in front 
of the star, receiving the stellar photons from behind, thought to happen during 
the superior conjunction of the compact object. In addition, the IC emission 
from secondary pairs or the occurrence of electromagnetic cascading, reprocessing 
energy above ~ 100 GeV down to lower energies, may affect the GeV variability as 
well, since the observable outcome of cascading is sensitive to the binary /observer 
geometry. 

Gamma-rays (TeV) 

The radiation processes contributing at these energies are the same as those 
contributing in the GeV range. Nevertheless, depending on the particle energy 
distribution and cross section behavior at high energies, some processes may be 
more or less significant in the GeV and TeV regimes. The nature of the emitting 
particles are also relevant. For instance, a synchrotron dominated particle energy 
distribution is much softer than a KN IC one, and electrons may be softer or 
harder than protons depending on the cooling conditions, although protons will 
in general yield a harder radiation spectrum at the highest energies for the same 
h. 

The variability at these energies can be produced by the same factors as in 
the GeV range. Regarding the angle-dependence of IC scattering, at TeV energies 
this process takes place in the KN regime. There is still some modulation of the 
emission, similar to some extent to (1 — cos(0)), but the deeper in the KN regime 
the IC scattering occurs, the smaller the modulation. At energies ^> TeV, there is 
basically no modulation. The changing angular dependence at different energies 
renders a complex spectral time behavior. There is additional modulation due to 
photon-photon absorption above the pair creation threshold, due to ~ (1— cos(0)) 
and electromagnetic cascading, and the pair creation threshold also depends on 
the photon-photon scattering angle (recall: eth = 2m^c^/(l — cos(0))eo). 



There could be additional IC components due to the presence of strong accretion disk and/or 
corona photon fields (e.g. Bosch- Ramon et al. 2006). 
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Fig. m gives an idea of the importance of the angle dependence of photon- 
photon absorption and IC scattering for the lightcurve. Unhke other sources of 
variabihty, this modulation occurs for sure if the star is providing the targets for 
these interaction processes. 

4. Final comments 

In this lecture, we try to show that to understand microquasars, first one can 
approach the problem using simple tools. This can be performed, partially, using 
the available phenomenological knowledge on the source properties (known ele- 
ments of the system/emitter and their relationships, observational information 
about the emission behavior). This allows the construction of phenomenological 
jet emission models which, at some point, fail to deepen in our understanding of 
microquasar non-thermal radiation. Therefore, in addition to this phenomeno- 
logical modeling, which allows us to set up the basic scenario and obtain some 
general idea of the global radiation behavior, a simple but systematic physical 
study source by source is also to be done. With this study, important constraints 
can be imposed on the particle acceleration efficiency, the impact of particle con- 
finement and propagation, and the plausible radiation and attenuation processes 
occurring in microquasars. 

Once the previous step is finished, it will be time to perform hydrodynami- 
cal and magnetohydrodynamical simulations, which are powerful heuristic tools 
to investigate further the complexity of the source structure (e.g. to understand 
the evolution of the jet, and its interaction with the stellar wind and large scale 
environment, etc.). Finally, using high quality data observations, provided by 
the next generation instruments, and the theoretical knowledge acquired with 
previous modeling and simulations, one can attack the problem of which mech- 
anism of particle acceleration is more suitable to explain the observed radiation, 
and the radiation mechanism itself, in the context of a more constraint physical 
scenario. 
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son regime, has the same dependence on the scattering angle as pair 
creation, and therefore it gets smaller around INFC and larger around 
SUPC. 



The physics of non-thermal radiation in inicroquasars 



19 



Albert, J. et al. 2007, ApJ, 665, L51 

Araudo, A. T., Romero, G. E., Bosch-Ramon, V., k Paredes, J. M. 2007, A&A, 
476, 1289 

Bednarek, W. k Giovannelli, F. 2007, A&A, 464, 437 
Benaglia, P. 2005, BAAA, 48, 127 

Biermann, P. L. & Strittmatter, P. A. 1987, ApJ, 322, 643 
Blumenthal, G. R. & Gould, R. J. 1970, RvMP, 42, 237 
Boettcher, M. k Schlickeiser, R. 1997, A&A, 325, 866 
Boettcher, M. k Dermer, C. D. 2005, ApJ, 634, L81 
Bogovalov, S. V. k Aharonian, F. A. 2000, MNRAS, 313, 504 
Bogovalov, S. V. k Kelner, S. R. 2005, ARep, 49, 57 

Bordas, P., Paredes, J. M., Bosch-Ramon, V. 2008, Int. J. Mod Phys. D, in 

press 

Bosch-Ramon, V. k Paredes, J. M. 2004, A&A, 417, 1075 
Bosch-Ramon, V., Romero, G. E., k Paredes, J. M. 2006, A&A, 447, 263 
Bosch-Ramon, V. 2007, Ap&SS, 309, 321 

Bosch-Ramon, V., Motch, C., Ribo, M., et al. 2007, AkA, 473, 545 
Bosch-Ramon, V., Khangulyan, D., k Aharonian, F. A. 2008, A&A, 482, 397 
Cheng K. S., Romero G. E. 2004, Cosmic Gamma-Ray Sources. Kluwer Aca- 
demic Publishers, Dordrecht 
Coppi, P. S. k Blandford, R. D. 1990, MNRAS, 245, 453 
Corbel, S., Fender, R. P., Tzioumis, A. K., et al. 2002, Science, 298, 196 
Donati, J.-F., Babel, J., k Harries, T. J. 2002, MNRAS, 33, 55 
Drury, L. O. 1983, Rep. Prog. Phys., 46, 973 

Fender, R. P., Belloni, T. M., k Gallo, E. 2004, MNRAS, 35, 1105 
Gould, R. J. k Schreder, G. P. 1967, PhRv, 155, 1404 
Hillas, A. M. 1984, ARA&A, 22, 425 

Kaufman Bernado, M. M., Romero, G. E., Mirabel, I. F. 2002, 385, LIO 
Kelner, S. R., Aharonian, F. A., k Bugayov, V. V. 2006, PhRvD, 74, 4018 
Khangulyan, D., Hnatic, S., Aharonian, F., k Bogovalov, S. 2007, MNRAS, 380, 
320 

Khangulyan, D., Aharonian, F., k Bosch-Ramon, V. 2008, MNRAS, 383, 467 
Mirabel I. F., Rodriguez L. F., Cordier B., Paul J., k Lebrun F. 1992, Nature, 
358, 215 

Li, Y., Wang, D. X., k Can, Z. M. 2008, 2008, A&A, 482, 1 
Mirabel, I. F., Dhawan, V., Chaty, S., et al. 1998, A&A, 330, L9 
Mirabel, I. F. k Rodriguez, L. F. 1999, ARA&A, 37, 409 

Orellana, M., Bordas, P., Bosch-Ramon, V., Romero, G. E., Paredes, J. M. 2007, 
A&A, 476, 9 

Owocki, S. P. k Cohen, D. H. 2006, ApJ, 648, 565 

Pacholkczyk, A. G. 1970, Radio Astrophysics (San Francisco: Freeman k Co.) 
Paredes, J.M., Marti, J., Ribo, M., Massi, M., 2000, Science, 288, 2340 



20 



V. Bosch-Ramon 



Perucho, M. & Bosch-Ramon, V. 2008, A&A, 482, 917 
Protheroe, R. J. & Stanev, T. 1987, Nature, 328, 136 

Reynoso, M. M., Christiansen, H. R., & Romero, G. E. 2008, Astrop. Phys., 28, 
565 

Ribo, M. 2005, ASP, 340, 421 

Romero, G. E., Torres, D. F., Kaufman Bernado, M. M., & Mirabel I. F., 2003, 
A&A, 410, LI 

Romero, G. E., Okazaki, A. T., Orellana, M., & Owocki, S. P. 2007, A&A, 474, 
15 

Rybicki, G. B. & Lightman, A. P. 1979, Radiative processes in astrophysics (New 

York: Wiley-Interscience) 
Semenov, V., Dyadechkin, S., k Punsly, B. 2004, Science, 305, 978 
Shakura, N. I. & Syunyaev, R. A. 1973, A&A, 24, 337 
Stelzer, B., Flaccomio, E., & Montmerle, T. 2005, ApJS, 160, 557 
Szostek, A., Zdziarski, A. A. 2007, MNRAS, 375, 793 
Usov, V. v., Melrose, D. B. 1992, ApJ, 395, 575 
Vila, G. S. & Romero, G. E. 2008, Int. J. Mod Phys. D, in press 



